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ABSTRACT 

Characterizing the physical properties of exoplanets, and understanding their formation and orbital 
evolution requires precise and accurate knowledge of their host stars. Accurately measuring stellar 
masses is particularly important because they likely influence planet occurrence and the architectures 
of planetary systems. Single main-sequence stars typically have masses estimated from evolutionary 
tracks, which generally provide accurate results due to their extensive empirical calibration. However, 
the validity of this method for subgiants and giants has been called into question by recent studies, 
with suggestions that the masses of these evolved stars could have been overestimated. We investigate 
these concerns using a sample of 59 benchmark evolved stars with model-independent masses (from 
binary systems or asteroseismology) obtained from the literature. We find very good agreement 
between these benchmark masses and the ones estimated using evolutionary tracks. The average 
fractional difference in the mass interval ~0.7 - 4.5 Mq is consistent with zero (-1.30 ± 2.42%), with 
no significant trends in the residuals relative to the input parameters. A good agreement between 
model-dependent and -independent radii (-4.81 ± 1.32%) and surface gravities (0.71 ± 0.51%) is also 
found. The consistency between independently determined ages for members of binary systems adds 
further support for the accuracy of the method employed to derive the stellar masses. Taken together, 
our results indicate that determination of masses of evolved stars using grids of evolutionary tracks is 
not significantly affected by systematic errors, and is thus valid for estimating the masses of isolated 
stars beyond the main sequence. 

Keywords: stars: fundamental parameters — stars: evolution — binaries: general — asteroseismology 


1. INTRODUCTION 

A complete understanding of the formation and evo¬ 
lution of planetary systems requires a detailed knowl¬ 
edge of the physical properties of the planets and their 
host stars. A subject of intense study within exoplane¬ 
tary science are the nature of the links between plan¬ 
ets and the stars they orbit. Planets form from the 
same molecular cloud material from which their host 
stars form, so planets can be thought of as the leftover 
relics of the star-formation epoch. Similarly, stars can be 
studied as relics of the environment from which planets 
formed long ago ( iLaughliii et alJl2004lFisc her fc Valentil 
1200, 5l [Kennedy fc KenvonI I2008I iJohnsonI I2008H . Rela- 
tionships between planetary systems observed today and 
their host star’s physical properties therefore provide im¬ 
portant constraints on theories of planet formation and 
orbital evolution. 

The two key physical properties that govern stellar 
evolution are mass and chemical composition, and these 
characteristics of host stars are thought to be linked to 
the formation and orbital evolution that produced the 
architectures of planetary systems seen today. Many in¬ 
dependent studies have shown that the occurrence of gi¬ 
ant planets around FGK dwarfs and post-main-sequence 
stars such as subgiant s strongly increases with the stellar 
iron abundance (e.g.. IGonzaledll997t ISantos et al.l[2004 

iFischer fc Valentill2005HJohnson et al.ll2010all . However, 

there exists much debate as to whether this correlation 
between metallicity and giant planet occurrence exists for 
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stars on the red giant branch (e.g, iJohnson et al.ll2010at 
Maldonado et alJ I20I3I IMortier et al.l 120131 IJofre et all 
2ni4IBeffert et aLll2ni 5ll or stars that host planets with 
masses less than that o f Neptune (e.g., Ghezzi et al.l 
l2010bl : iSousa et al.l 120111: iBuchhave et al.l 201411 . Evi- 
dence for possible connections between the presence of 
planets and other peculiarities in the chemical abun- 
danc es of their stellar hosts are still ambiguous (see , 
e.g.. lAdi bekyan et al.l 120121 12014 iFigiieira et al.l 12014 
iTeske et al.ll20l4 and references therein). 

Also still unclear is the role of stellar mass on the 
formation and evolution of extrasolar planets. Planet- 
search programs based on the radial velocities (RVs) 
technique have mainly targeted Solar-type stars be¬ 
cause they are relatively numerous and bright in the 
Solar neighborhood, and their spectra have a consid¬ 
erable number of unblended, narrow metallic lines, al¬ 
lowing more precise measurements of both the Doppler 
shifts and stellar pr o perties (e.g. iGalland et al.l 120051 : 
iLagrange et al.l 120091 : iBecker et al.l I2015IL One direct 
consequence of this selection criterion is a relatively 
small range of stellar masses included among the tar¬ 
gets of the largest existing planet-search programs (~0.7 
-1.4 Mq). The extension of this interval to lower-mass M 
dwarf s was a natural extens i on of the exop l anet searches 
(e.g., IJohnson et al.l l2007al: iBonfils et al.l I2013H . given 
these stars are more numerous in the Galaxy and al¬ 
low robust detections of Earth-like pl anets with the cur¬ 
rent instrumental ca pabilities le.g.. iMavor et ahl I20l4 
iQuintana et al.ir2014ll . Nevertheless, the monitoring and 
characterization of M dwarfs has been challenging be¬ 
cause they are faint and their optical spectra are com- 
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plete l y dominated by m olecular bands (e.g. iBean et’aLl 
l2006t iManess et ciI1l2007D . Dedicated ongoing and future 
suryeys will progressiyely oyercome thes e issues (e.g., 
iBerta et al.ll2013t iQuirrenbach et ani2014li . 

A completely different difficulty is encountered for 
planet searches around massive stars on the main se¬ 
quence (>1.5 Mq). The hotter effective tempera¬ 
tures and higher rotational velocities of early F and A 
stars produce a spectrum with fewer and much broader 
lines, hindering the measurement of precise radial ve- 
locities and the reliable detection of ex trasolar planets 
(|Galland et al.l [2?i?)5 iBecker et al.l 120151 1. However, as 
these massive stars evolve off the main sequence towards 
the red giant branch (RGB ), they cool and slow down 
(|do Nascimento et al.l 1200011 . As a result, their spec¬ 
tra then show many narrow lines which allow the mea¬ 
surement of Doppler shifts caused by extrasolar planets. 
Therefore, analyzing the evolved counterparts of main- 
sequence F- and A-type dwarfs is currently the best al¬ 
ternative to investigate planet occurrence around more 
massive starfl 

To this end, more than 1000 e volved stars have been 
monitored by several surveys Isee iNiedzielski et al.ll2015l 
and references therein), resulting in the detection o f more 
than 100 planets around them (iJofre et al.l I20l4) . Al¬ 
though this sample is relatively small when compared to 
the entire sample of discovered planets, some interesting 
properties already started to emerge, such as a paucity 
of plane ts at short orbital distances and large eccen¬ 
tricit ies (jSato et al.ll2007t[Bowler et al.ll20ir)l : lJones et al.l 
120141 and t he possible lack of a giant planet - metallicity 
correlation (iGhezzi et al.ll2010al: iMaldonado et al.ll2011^ 

IJofre et al.ll20l4 but see also iReffert et al.ll2015ll . 

Of particular interest to this work is the higher oc¬ 
currence rate of J ovian planets a round more massive 
stars measured by iJohnson et al.l (j2010al l. The analy¬ 
sis of 1194 stars with masses in the range 0.2 - 2.0 Mq 
not only provided additional confirmation to the well- 
established giant planet-metallicity correlation, but also 
revealed that the occurrence rate of these planets in¬ 
creases linearly from ^3% for M dwarfs to ~14% for 
A stars. This is a very important result because it adds 
new constraints to the planet formation theories and also 
provides guidance to ongoing and future surveys moni¬ 
toring evolved stars (rad ial velocities and transit s) or A 
dwarfs (direct imaging - iCreop fc Johnsonll20lH) . 

The determination of masses for these isolated evolved 


stars relied on the comparison of observed properties 
with stella r evolutionary tracks ([Johnson et al.l l2010al 
|2011L[20T^ . The different models usually adopted in this 
method were subjected to numerous improvements and 
validation tests over the last few decades, leading to more 
precise predictions of the stellar parameters and also a 


lution of stars (e.g., 

Vandenberd 119851: lAndersen et al. 

1988 

; lAndersenI II991 

: Andersen et al.l 11991: IPols et al. 

1997 

; Demaroue et al. 2004 Torres et al.l 20101: 

Bressan et al.l 2012 

: Brogaard et al.l 120121: IGarcia et al. 

2014 ITorres et al. 

20151). For field stars on the main 


sequence, this technique for determi ning masses gener- 
ally provides accurate results (e.g., iTorres et al.l 1201^ 
that are a lso fairly consistent with other independent 
estimates (jPinheiro et al.l 12014 but note the possible 
issues with spectroscopic masses for stars with M > 1.2 
Mq). 

The reliability of the application of this method to sub¬ 
giants and giants is, however, more uncertain. Although 
stellar evolution models were able to successfully describe 
some binary systems with a t leas t one e volved com¬ 
pone nt (e.g., lAndersen et al.l 119881 119911 : iTorres et al.l 
l2015f l. there were suggestions that the masses of planet¬ 
hosting subgiants and giants could have been overesti¬ 
mated (by up to 50%) due to systematic errors o n their 
atmospheric parameters or the models themselves (iLlovdl 
1201 IL [2(TT^ iSchlanfman fc WinnI 1201311 . Additional evi¬ 
dence of possibly overestimated masses for evolved stars 
wi th planets (by up to 100%) were recently presented 
bv iSousa et al.l 1)20151 ) (but note that a possible explana¬ 
tion f or some of the problemat ic stars is given in Section 
4.2 of iTakeda fc Taiitsul 1201511 . In spite of the growing 
number of evolved stars with precisely determined masses 
(through orbital solutions in binary systems or asteroseis- 
mology of field or cluster stars), the evolutionary tracks 
continue to be main method to determine this funda¬ 
mental parameter for single subgiants and giants. Thus, 
possible errors in the results would have implications on 
many different studies, ranging from the formation and 
architectures of planetary systems to Galactic chemical 
and dynamical evolution, and have to be carefully inves¬ 
tigated. 

We and others have been conducting a study to pre¬ 
cisely constrain the masses of evolved stars using differ¬ 
ent input parameters and techniques. The first results of 
this effort for the b right nearby subgiant star HD 185351 
were presented bv IJohnson et all (|2014D . In this work, 
we check the accuracy of the evolutionary tracks method 
using a sample of benchmark subgiants and giants with 
accurate masses. This sample and the literature masses 
for their masses are described in Section [2l while the de¬ 
termination of corresponding stellar masses from evolu¬ 
tionary tracks is discussed in Section |3l In Section IH we 
compare the two sets of masses and show that the evolu¬ 
tionary tracks method do not seem to have any system¬ 
atic errors that would overestimate the masses of evolved 
stars. Finally, our concluding remarks are presented in 
Section m 


2. SAMPLE OF BENCHMARK STARS 

Our sample of benchmark stars consists of 59 subgiants 
and giants with precise masses determined with meth¬ 
ods that are based on minimal assumptions and physical 
modeling. In this study, we focus on dynamical masses 
measured for stars in binary systems and asteroseismic 
masses calculated with scaling relations and the direct 
method (i.e., without the usage of [Fe/H] and grids of 
models t o further constra in the evolutionary parameters; 
see, e.g.. lGai et aLll^llll . 


^ Direct imaging is a promising technique to detect giant planets 
around B- A stars, but so far only a handful of them were found 
(see, e.g., ICrePD fc .Tohnsor] 120111 : (Nielsen et al.l 120131 and refer¬ 
ences therein) 


2.1. Stars in Binary Systems 

We performed an extensive literature search for de¬ 
tached binary systems in which at least one of the com¬ 
ponents is an evolved star. Most of them were found with 
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the help of DEBCal[3 (|Southworthll2014[l . Only subgiants 
and giants that have available values for effective tem¬ 
perature (Teff), inetallicity ([Fe/H]), V magnitude and 
parallax (tt) or distance (d) were selected. These are 
the input parameters necessary to obtain masses with 
the evolutionary track method (see Section [3]). We re¬ 
moved any stars with metallicities determined through 
the best fit to evolutiona. r y tracks or isochrone s (e.g., 
iSandberg Lacv et al.ll2012l: iRataiczak et al.l[2013ll in or¬ 
der to avoid dual dependencies on the model tracks, e.g. 
using [Fe/H] derived from one evolution model grid to 
interpolate onto another model grid. 

Our selection criteria resulted in a sample of 26 evolved 
stars in 16 binary systems in the Milky Way (MW), 
Small Magellanic Cloud (SMC) and Large Magellanic 
Cloud (LMC). These stars are shown in Tables [1] and 
[3] along with their fundamental parameters, respective 
uncertainties and literature references. In some cases, 
no uncertainty was available for the V magnitudes and 
we adopted cr(V) = 0.01 for MW stars because this 
is the typical accuracy of the Hipparcos measurements 
IPerrvman et al.l (|1997ll . For LMC stars, we adopted a 
conservative uncertainty cr{V) = 0.02 because the accu¬ 
racy of photometric ca librations in the OG LF project is 
better than 0.02 mag l|Udalski et al.l[2058ll . We also as¬ 
sumed E{B — V) = 0.000 for all stars closer than ^^60 pc. 
Their positions in a HR diagram can be seen in Figure 
[TJ Dynamical stellar masses and radii were directly de¬ 
termined in the original source papers with precisions of 
<3% and <7%, respectively, through the analysis of ra¬ 
dial velocities (RVs) and light curves. The masses range 
from ^1.2 to ^4.5 Mq. 

2.2. Stars with Asteroseismic Parameters 

We conducted a separate literature search for evolved 
single stars with measurements of the two global astero¬ 
seismic parameters: the frequency of maximum power 
r'max and the large frequency separation Ah'. Once 
again, we kept only those stars which also had reported, 
evolution-model-independent values for Teg, [Fe/H], V 
mag nitude and parallax (tt) or distance {d) (see Section 
[2T]) . The last requirement, in particular, dramatically 
reduces the number of stars suitable for our study. 

We calculated asteroseismic masses and radii with 
the direct method, using published values of 
Au and T^g as input t o the scaling relations (e.g., 
iKieldsen fc Beddinglll995ll written as: 




Figure 1. HR diagrams showing the positions of stars in our 
sample. The upper panel shows the MW binaries (black points) 
and field stars with asteroseismology (red points). The lower panel 
presents the eclipsing binaries (EBs) from LMC (blue points) and 
SMC (green points). In all panels, the solid lines represent PAR¬ 
SEC evolutionary tracks l|Bressan et al.ll2Q12l ^ for 1.0 Mq, 1.5 Mq, 
2.0 Mq, 3.0 Mq and 4.0 Mq, and representative metallicities for 
the depicted samples. 


M / \ Y Teff 

Mq-[ 

^max,0 

Rq yiZjjjax,© J \yAh'Q J \yTeg^0 J 

The reference solar asteroseismic paramete rs = 

3090 ± 30 ^Hz and Aiz© = 135.1 ± 0.1 /iHz (jHuber et al.l 
[Ml were adopted in the above equations, as well as 
the canonical value T^g^Q = 5777 K. We note that, while 
these results are independent from stellar evolution mod¬ 
els, they do rely on scaling of solar p-modes to stars of dif- 

^ http://www.astro.keele.ac.uk/jkt/debcat/ 


ferent masses, metallicities and evolutionary states. Re¬ 
cent studies have shown that the usage of the scaling re¬ 
lations for inetal-poor or evo l ved stars require som e cau¬ 
tion (e.g., jHIilioIer^lloil; [Fpstein^E^llinU- How¬ 
ever, many different tests indicate that these scaling re¬ 
lations provide masses and radii with accuracies ~10% 
and ^5%, respectiv ely (see discussion in Section 2.4.1 of 
[Johnson et ahllMll . No corrections were applied to the 
scaling relations because there is still no consensus which 
of th e proposed ones works best (e.g., iBrogaard et al.l 
(Ml . Nevertheless, the shifts in M and R that would 
be caused by adopting the available corrections are dis¬ 
cussed in Section m 

We should also note that Equations [T] and [2] con¬ 
tain the stellar effective temperature, which in turn can 
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be determined using model-dependent techniques (for 
instance, spectroscopic Tefj depends on model atmo¬ 
spheres). However, we note that the dependence of the 
stellar radius on in Equation[2]is weak relative to the 
other terms. Moreover, the uncertainties on the values of 
masses and radii caused by typical errors on the effective 
temperatures (^80 K or ~2% considering a typical Teff= 
^4900 K) are negligible (of order 1%). 

As one of our goals is to define a statistically significant 
sample of benchmark stars, we decided to remove stars 
for which the uncertainties in the asteroseismic masses 
were higher than 20%. After this cut, we obtained a 
subsample of 33 evolved stars with asteroseismic param¬ 
eters and masses with typical precisions ranging from 5% 
to 20%. These relatively large errors could be decreased 
if we used a grid-based asteroseismic method (jGai et al.l 
unni), in which the stellar properties predicted by a grid 
of stellar evolution models are converted to the astero¬ 
seismic parameters z^max and Aiy, and then compared to 
the observed values. But the goal of this work is to avoid 
any model dependencies in the reference parameters of 
the benchmark stars sample. 

Our asteroseismic benchmark stars are shown in Tables 
[Hand [21 As for the binaries, u(V) =0.01 was adopted for 
stars with no uncertainties available for V. We assumed 
E{B — V) = 0.000 for all stars closer than ^60 pc as well 
as 10 eight more distant stars for which no values were 
available in the literature. When no uncertainties were 
not reported for the asteroseismic parameters, arbitrary 
errors of 5% and 2% were adopted for i/max and Ai/, 
respectively ('following iBruntt et al.ll2010ll . The positions 
of the stars in a HR diagram is shown in Figure [1] and 
their masses range from roughly 0.7 Mq to 4.0 Mq. One 
star (KIC 8410637) is common to both the binary and 
asteroseismic samples. In this case, we adopt the more 
precise dynamical stellar mass and radiufl 


2.3. Final Sample 

Our final sample consists of 59 benchmark stars with 
precise dynamical or asteroseismic masses. According 
to model-dependent classifications provided in the ref¬ 
erence papers or to Figure [H we can see that most of 
our stars are on the red giant branch (RGB) or on the 
red clump (RC), which are also the regions on the H-R 


cus on evolved stars 

e.g.. 

Hatzes et al.l[2005[: [Sato et al. 

2005 

Lovis & Mavod 

2007 

: Lee et al.l 20111: Omiva et al. 

2012 

NiedzieLski et al.l[20^ 

51 [R.effert et ali[2015D. All of 


them have the required parameters to allow an indepen¬ 
dent determination of their masses using stellar evolution 
models (see section El). 

As explained in the previous sections, some stars were 
removed from our sample because they did not fulfill one 
or more of our selection criteria. For completeness, those 
stars are listed in Table E] along with the selected refer¬ 
ences and the reasons for their exclusion. It is also worth 
noting that we may have overlooked a few stars despite 
our best efforts to conduct a thorough literature searcl: 0 . 
However, we consider that our current sample is statis¬ 
tically robust, and the addition of a few objects will not 


likely affect our conclusions. 

3. STELLAR MASSES FROM EVOLUTIONARY TRACKS 

We derived model-dependent stellar masses for all 59 
of our benchmark stars using stellar evolution model 
tracks. In this method, measured stellar properties 
(typically, Tes, log-luminosity log^Q (L/Lq), and metal- 
licity [Fe/H]) are compared with the properties pre¬ 
dicted by evolution models calculated for different ini¬ 
tial masses and chemical compositions (e.g.lGirardi et al.l 
120021: [Potter et al.ll2008HDemarane et al.l^OSl ). A orob- 

abilistic analysis then provides the best match between 
the observed and theoretical properties of the star, al¬ 
lowing the determination of its mass (as well as radius 
and age). This method is known to work well for solar - 
type stars on the main sequence (e.g., [Torres et al.ll20ldl) . 
but has not been extensively tested on f or stars in ad¬ 
vanced evo lutionary states (however, see [Johnson et al.l 
IMjLIMI . Testing these model grids using a large sam¬ 
ple of evolved stars is the primary goal of this study. 

We adopted the implementation of the method con¬ 
tained in the PARAM code, which was kindly provided 
by Leo Girard0. Estimates of the mass, radius, log 5 and 
age are obtained through a Bayesian estimation method 
using different priors , input parameters and isochrones 
l|da Silva et al.[[200^ . We adopted the default options 
for the Bayesia n priors and s e lected the new PARSEC 
isochrones from [Bressan et Ml ([20I2| ). We also chose the 
analysis of stars with known parallaxes, for which we pro¬ 
vided the additional input parameters Teff, [Fe/H] and V 
magnitude (see Table[H). Extinction corrections Ay were 
applied to the V magnitudes of all stars farther than 60 
pc for which a reddening value E{B — V) was available in 
the literature. The standard relation Ay = 3AE{B — V) 
was adopted for the conversion of reddening to H-band 
extinction. 

It is worth noting that the Padova models (in their 
most recent version, PARSEC: [Bressan et al.l[2Mll were 
chosen because its wide usage in the literature, includ- 
ing the estimation the prop erties of evolved star s (e.g. , 
Johnson et al.[[2?)^ [2007b l: [Pollinger et al.[ [20071 [20091 
de Medeiros et al.[ [20091 [Ghezzi et al.l [20I0aU bh. Tests 
with other grids of evolutionary tracks are beyond the 
scope of the present work, but will be adressed in a fu¬ 
ture contribution using the same sample of benchmark 
stars. 

The masses derived with the evolutionary tracks are 
shown in Table [2| along with the other output parame¬ 
ters of PARAM and their respective uncertainties. For 
the giant star ylnd, the error on [Fe/H] was arbitrar¬ 
ily increased to 0.20 dex because the PARAM code 
was not returning valid solutions with the original value 
(cfFe/H] = 0.07). After some tests, it was clear that 
PARAM was having issues due to the relatively low 
metallicity and small original uncertainty. The preci¬ 
sion of the masses obtained with the evolutionary tracks 
method varies from ^2% to ~29%, with a typical (me¬ 
dian) value of ~I0%. The relative uncertainties on the 
masses are similar for the samples of binaries and astero¬ 
seismic targets and increase for stars with lower effective 


^ The dynamical and asteroseismic masses, radii and surface 
gravities agree within 1.6 (t (see Table [2j. 

^ Please, contact the authors if you have identified such a case. 


° For the web interface maintained by Leo Girardi at the Osser- 
vatorio Astronomico di Padova, visit http://stev.oapd.inaf.it/cgi- 
bin/param. 
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temperatures and larger errors on their metallicities. 

4. RESULTS AND DISCUSSION 
4.1. Masses 

Our comparison between model-independent masses 
and the corresponding values derived from evolutionary 
tracks (see Section [3|) is shown in Figure [2j We can see 
there is an overall good agreement and no systematic 
offsets for the entire mass range (~0.7 - 4.5 Mq), even 
though our results are based on a heterogeneous data set. 
The average absolute and percentage differences between 
evolutionary track and reference masses for the entire 
sample and the two subsamples (binaries and asteroseis- 
mic targets) are shown in Table H] We can see that the 
global average differences are consistent with zero within 
the uncertainties, with 40 and 33 stars (i.e., 68% and 
56% of the sample, respectively) showing agreements be¬ 
tween the two mass estimates within 20% and 10%, re¬ 
spectively. The separate results for the binary and as- 
teroseismic samples are somewhat different. While the 
average mass difference for the former is zero within the 
errors, it reveals that the evolutionary track results are 
slightly underestimated for the latter. We therefore see 
no evidence that the models overpredict the masse s of in¬ 
dividu al st ars, in contrast to the concern s raised bv iLlovJ 
(|2013ll and iSchlaufman &: Wind (|2013f) . Note, however, 
that the stars analyzed here and in those works are in 
different evolutionary stages. 

We performed weighted linear fits to the fractional 
residuals of the whole sample as well as the asteroseis- 
mic and binary subsamples individually. The weights 
were calculated as where atotai is the fractional 

uncertainty obtained with the propagation of errors for 
the difference of the two mass estimates. The results 
of the fits are listed in Table |4] and demonstrate that a 
slight trend is observed for the asteroseismic subsample, 
with the model grids underpredicting stellar masses by 
an amount that increases with mass. This is also clear 
from the lower panel in Figure [2] 

The results for the asteroseismic targets deserve some 
further discussion. Thirteen out of 33 stars have frac¬ 
tional differences larger than 20%: Arcturus, /3Aql, 
S Eri, rj Ser, HD 170008, HD 170231, HD 178484, 
KIC 4044238, KIC 6101376, KIC 7909976, K I C 8508 931, 
KIC 8813946, and ^ Hya. iKallinger et al.l (j201fl[ l de¬ 
rive a mass M* = 0.80 ± 0.20 Mq for Arcturus (using 
Ai/ and an interferometric radius), which agrees better 
with the value derived here from the evolutionary tracks 
(M* = 0.98±0.06 Mq ). The adoption of the larger aster¬ 
oseismic mass would reduce t h e frac tional difference from 
45% to 22%. iBruntt et al.l (j2010fl provide alternative 
masses for /3 Aql, <5 Eri and rj Ser: M* = 1.26 ±0.18 Mq 
(calculated from iZmax, Teff and interferometric radius), 
M* = 1.33 ±0.07 Mq (calculated from Av and interfer¬ 
ometric radius) and M* = 1.45 ± 0.21 Mq (calculated 
from i^max, Teff and luminosity), respectively. All these 
literature masses are in better agreement with the ones 
obtained from evolutionary tracks and their usage would 
decrease the fractional differences from 31%, 36% and 
30% to 10%, 12% and 15%, respectively. 

KIC 8508931, has an alternative determination of the 
globa l asteroseismi c para meters. Adopting the results 
from iHekker et al.l (|2011h (but using the same Tgff as 


before and assuming uncertainties of 5% in t'max and 
2% in Aiz), we obtain a mass M* = 2.54 ± 0.43 Mq, 
which reduces the discrepancy fr om 26% to 12%. For 
KIC 8813946, iHuber et al.l (|2012ll found a conflict be¬ 
tween the radii calculated for this star using different 
constraints (asteroseismic parameters and angular diam¬ 
eter from interferometry coupled with the parallax) and 
argue that this discrepancy could be indica tive of a prob¬ 
lem w ith the revised Hipparcos parallax (|van Leeuwenl 
I2007D . This problem would, in turn, affect our mass de¬ 
termined from the evolutionary tracks. Finally, the star 
^ Hy a also has an alterna tive mass M* = 2.89 ±0.23 Mq 
from IBruntt et al.l (I2010D . calculated from t'max, Teff and 
its interferometric radius. The adoption of this mass 
would improve the agreement, with a decrease in the 
fractional difference from 23% to 16%. Unfortunately, 
we were not able to find alternative parameters for the 
stars HD 170008, HD 170231, HD 178484, KIC 4044238, 
KIC 6101376, and KIC 7909976. 

If we replace the masses of the problematic stars with 
the alternative values mentioned above for six of them, 
the average differences become: (AM) = -0.22 ± 0.07 
Mq and (AM/M/je/.) = -7.18 ± 2.80%. We can see that 
the fractional average difference became more negative 
because the large positive values of Arcturus, /3 Aql and 
5 Eri were decreased in the first case and replaced by 
negative offsets for the other two stars. The absolute 
average difference did not change significantly and no 
improvement in the standard deviations of the mean was 
observed. The parameters of the linear fit, on the other 
hand, were slightly improved: A = —7.98 ± 2.34, B = 
2.76 ± 5.07, = 0.274 and a = 0.77. Therefore, it is 

clear that part of the trend observed for the asteroseismic 
targets was caused affected by these six stars with large 
discrepancies and, with the exception of rj Ser, located at 
the edges of the mass interval covered by the sample. The 
remaining trend in the residuals can be mostly attributed 
to the other seven stars with large deviations and no 
alternative parameters to be tested. 

Corrections to the s calin g relations have been proposed 
bv I White et al.l l)2nilh and iMosser et al.l (120131) . but fur¬ 
ther tests are still necessary to confirm their accuracy. 
As a test, we applied the latter to the 33 stars in our 
asteroseismic sample (the former is valid only down to 
Teff = 4700 K and some stars in the sample are cooler 
than this limit) and observed a typical decrease in the as¬ 
teroseismic masses of ^5%, which increases the average 
fractional difference {AM/Mnef.) from -4.29% to 0.47%. 
Although the agreement is better, we decided to keep our 
original masses because there is still no consensus in the 
literature regarding the validity of the proposed correc¬ 
tions to the scaling relati ons or which one works best 
(e.g., iBrogaard et al.l[2014[) . 

The exercise of testing alternative parameters could 
not be performed for the six binaries that have percent¬ 
age differences larger than 20%: RT CrB B, OGLE SMC 
130.5 4296 A, OGLE LMC-ECL-01866 B, OGLE LMC- 
ECL-03160 A, OGLE LMC-ECL-09660 B and OGLE 
LMC ECL15620 A. To the best of our knowledge, there 
are no alternative reference or input parameters that we 
could use to check if a better agreement is obtained. Also, 
we were not able to observe any clear pattern in these 
stars with larger deviations. They are all from differ¬ 
ent binary systems and their companions seem to yield 
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Figure 2. Comparison between evolutionary track and reference masses. Milky Way binaries, LMC EBs, SMC EBs and field stars 
with asteroseismology are shown by black, blue, green and red points, respectively. Black solid lines represent a perfect agreement. The 
upper panel shows a direct comparison of the masses. The lower panel shows the fractional difference between the mass estimated from 
the evolutionary tracks and the reference mass. We can see there is a general good agreement between the masses and no significant 
overestimate of the values by the evolutionary track method. 


at least reasonable results (percentage differences lower 
than 15%). Moreover, they are not all the primary or 
secondary components in their systems. Thus, it is not 
feasible to attribute any problems to the analysis or the 
parameter determination for specific systems. The only 
detail worth noting is that four of the six stars are from 
the LMC, i.e., are extragalactic in origin. 

We also investigated if the larger fractional differences 
could be caused by the degeneracy between the two pos¬ 
sible evolutionary stages for most stars in our sample: 
red giant branch (RGB) or the red clump (RC). This de¬ 
generacy arises from the identical absolute magnitudes 
(luminosities) along the ascent of the RGB and stars on 
the horizontal branch, or red clump region of the H-R 
diagram. Uncertainties in the effective temperatures and 
metallicities of these high-luminosity giants therefore re¬ 
sults in a do uble-peaked probabil ity distribution function 
(PDF) (e.g.. lda Silva et aT1l2006ll and this pattern can be 
indeed observed for many of our stars. A possible con¬ 
sequence of this issue could be the choice of an incorrect 
mass associated with just one of the peaks (for instance, 
the mass corresponding to the maximum of the PDF), 


creating a systematic offset for at least some of the stars 
in the sample. 

Our results do not show such offsets: evolutionary 
tracks are able to correctly recover the reference masses 
with an average accuracy better than 10% and no sig¬ 
nificant systematic trends. This good agreement is re¬ 
lated to the way PARAM determines the best solution: 
it calculates mean values for the parameters from their 
PDFs, instead of choosing the most likely values, for ex¬ 
ample. Although this increases the uncertainties on the 
results, it avoids attributing the parameters to a star 
given a choice of a specific evolutionary state. There¬ 
fore, the possible degeneracies are taken into account in 
our analysis and do not to seem to be the reason for 
the most discrepant cases. We should also note that the 
differences between the two peaks of most of the double- 
peaked PDFs are within 68.3% confidence interval of the 
mean value. 

As a final check on our results, we checked if the mass 
differences were a function of any of the parameters that 
were used as input to the evolutionary tracks method 
(Teff, [Fe/H], V magnitude and parallax). As can be 
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seen in Figure [31 no systematic trends are found for any 
of the parameters. The higher correlation coefficients 
of the weighted linear fits are 0.020 for the entire sample 
(in Teff), 0.014 for the binaries (in [Fe/H]) and 0.124 (in 
parallax) for the asteroseismic targets. 

4.2. Radii 

Our analysis of stellar radii closely follows that of the 
previous section. The comparison between reference radii 
and the corresponding values estimated from evolution¬ 
ary tracks (see Section |3|) is depicted in Figure [H Al¬ 
though the overall agreement is good, we can see a global 
small offset (^5%, see Table II, with the evolutionary 
tracks underestimating the radii of most stars. There is 
also a clear difference between the samples of binaries 
and asteroseismic stars. While the latter have a larger 
and mostly constant offset with a significant dispersion 
around it, the former have a smaller dispersion around 
a small offset that slightly increases with increasing ra¬ 
dius. It is worth noting though that these offsets are 
well within the uncertainties obtained in the radii from 
evolutionary tracks, as is also clear from Figure |H 

Besides these interesting features, once again we can 
find asteroseismic targets with percentage differences in 
radius larger than 20%: ^ Aql, HD 169730, HD 170174, 
HD 170231, KIC 8508931 and KIC 8813946. Four of 
them also had problems in their masses. For /3 Aql, our 
evolutionary track radius (3.25 ± 0.13 Rq) is in very 
go od agreement with th e interferometric val ue presented 
bv iBruntt et al.l (|2010ll (3.21 ± 0.13 Rq). iHuber et al.l 
(j2012tl provide an interferometric radius for KIC 8813946 
(R = 12.0 ± 1.2 Rq) which agrees very well with the 
value derived from the evolutionary tracks (12.0 ± 1.3 
Rq). This could suggest a problem in the asteroseis¬ 
mic global parameters for this star, but we should re¬ 
call that the parallax is used as an input for both the 
interferometric and evolutionary track radii. Thus, we 
can not d iscard a possible is sue on the parallax (as sug¬ 
gested bv IHuber et al.l[?012[l . which could be leading to 
two agreeing, but incorrect radii. Finally, adopting for 
KIC 85089 31 the alterna t ive gl obal asteroseismic param¬ 
eters from iHekker et al] (|2011ll (but using the same 
as before and arbitrary uncertainties of 5% in i/max and 
2 % in Ai^) has a small impact on the radius, reducing the 
percentage difference from 23% to 19%. Unfortunately, 
we were not able to find alternative parameters for the 
stars HD 169730, HD 170174 and HD 17023T _ 

Cor rections to the scaling relations from iMosser et akl 
(j2013H were again applied as a test to asteroseismic sam¬ 
ple. We observed a typical decrease in the asteroseismic 
radii of ^2.5%, which increases the average fractional 
difference (AR/Rnef.) for the asteroseismic targets from 
-6.58% to -4.31%. and improves the agreement between 
the two sets of radii. Despite this improvement, we de¬ 
cided to keep our original radii for the same reason as 
explained for the masses. 

As in Section l4Tl we investigated if the offsets between 
the two sets of radii were correlated with Tgff, [Fe/H], V 
magnitude and parallax. No systematic behaviors were 
found in any case, with maximum correlation coefficients 
R? for the weighted linear fits equal to 0.044 for the entire 
sample (in Tes), 0.204 for the binaries (in [Fe/H]) and 
0.133 (in parallax) for the asteroseismic targets. 


4.3. Surface Gravities 

The last of the three parameters that was derived from 
both model-independent and -dependent methods is the 
surface gravity. The comparison between the two sets 
of results is shown in Figure |5| The agreement is re¬ 
markable, with average absolute and percentage differ¬ 
ences 0.01 ± 0.01 dex and 0.71 ± 0.51%. Such a good 
agreement is also observed for the subsamples of bina¬ 
ries and asteroseismic targets. The more robust results 
for the surface gravity are expected since this parameter 
(rather than the mass or radius) is directly connected 
to the astero seismic observables through the scaling re¬ 
lations le g.. iKie ldsen fc Beddind 119951 1. Moreover, as 
noted bv iGai et al.l ([2011 1. the errors on asteroseismic 
masses and radii have a strong positive correlation that 
compensate each other to produce small uncertainties in 
log g. For this reason, the usage of masses and radii de¬ 
rived from cor rected scaling relati ons (following the pre¬ 
scription from iMosser et al1l2013fl to calculate corrected 
surface gravities leads to virtually equal values of log g. 

As was the case for both mass and radius, no system¬ 
atic trends were observed in the percentage differences 
between the two sets of surface gravities as a function of 
any of the input parameters (Teff, [Fe/H], V magnitude 
and parallax). The higher correlation coefficients R^ for 
the weighted linear fits were 0.039 for the entire sample 
(in [Fe/H]), 0.041 for the binaries (in [Fe/H]) and 0.107 
(in Teff) for the asteroseismic targets. 

4.4. Ages 

The ages for all our benchmark stars were also given 
as an output by PARAM (see Table [2|), but we were not 
able to compare them with the literature values because 
these were also determined using theoretical isochrones. 
Although not all ages were helpful in assessing the per¬ 
formance of the evolutionary track method, the ones de¬ 
rived for binary systems which have both members in 
our sample can provide important information about the 
consistency of our results. The determination of ages 
for a binary system usually requires that both stars are 
described by a single isochrone at a given age, which is 
a reasonable assumption considering the stars were pre¬ 
sumably born at the same time from the same molecular 
cloud. However, we did not impose this constraint here 
and analyzed members of a given system as if they were 
individual isolated stars. Thus, a good agreement be¬ 
tween the ages derived for the stars in a binary system 
could be regarded as a consistency check for the results 
obtained with the evolutionary track method. 

Our sample contains 10 systems for which both stars 
were analyzed: 7 from the LMC, 2 from the SMC and 
one from the MW. The comparison between the ages 
of the members is shown in Figure [6l We can see the 
overall agreement is good, with only three systems not 
having ages that agree within la: OGLE SMC 130.5 
4296, OGLE LMC-ECL-09660 and OGLE LMC-ECL- 
26122. For OGLE SMC 130.5 4296, the A component has 
a mo del-dependent mass 37% higher than the dynamical 
one, which is consistent with its lower age relative to the 
B component. A similar situation is observed for OGLE 
LMC-ECL-09660, with the only difference that it is the 
B component that has a larger model-dependent mass 
(by 27%). For the system OGLE LMG-ECL-26122, both 
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Figure 3. Percentage difference between evolutionary track and reference masses as a function of Tgff (upper left panel), [Fe/H] (upper 
right panel), V magnitude (lower left panel) and tt (lower right panel). Symbols and lines are the same as in Figure [2] We observe no 
systematic trends for any of the parameters. 


components present a good agreement between the two 
sets of mass (5% for A and -4% for B). However, the 
small offsets in opposite directions are enough to cause 
the age discrepancy, especially because of the relatively 
small errors. 

The average absolute and fractional differences be¬ 
tween the ages of the components (in the sense A-B) are 
-0.15 ± 0.13 Gyr and -15.08 ± 13.80 % (where the errors 
are standard deviations of the mean). A weighted linear 
fit to the residuals does not show any systematic trends 
(i?^ = 0.033). Thus, the consistency check with the ages 
provides additional support to the good performance of 
the evolutionary track method for evolved stars. 

5. CONCLUSIONS 

The determination of masses for evolved stars (sub¬ 
giants and giants) is still a matter of debate due to recent 
claims in the literature that the interpolation of observed 
properties (such as Tgg, log(L/L 0 ) and [Fe/H]) in a grid 


of evolutionary tracks could yield systematically over¬ 
estimated results. The accuracy of this method has not 
been extensively tested in this region of the Hertzsprung- 
Russell diagram. In order to address this issue, we com¬ 
piled a sample of 59 benchmark stars from the litera¬ 
ture, which includes 26 members of binary systems (in 
the MW, SMC and LMC) and 33 objects with global 
asteroseismic parameters. 

We determined model-dependent masses for all stars 
using the PARAM model-grid interpolation code and 
then compared with model-independent values coming 
from dynamical (binaries) or asteroseismic (single stars) 
mass measurements. We observed a very good agree¬ 
ment for the entire mass interval between '^0.7 - 4.5 
Mq, even though a heterogeneous data set was used in 
the study. The average absolute and fractional differ¬ 
ences between model-dependent and reference masses are 
—0.10 ±0.05 Mq and —1.30 ±2.42%, respectively (where 
the errors are standard deviations of the means). No sig- 
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Figure 4. Comparison between evolutionary track and reference 
radii. Symbols, lines and panels are the same as in Figure [2] We 
can see there is a general good agreement between the radii with 
small offsets between the two sets of measurements. 
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Figure 5. Comparison between evolutionary track and reference 
surface gravities. Symbols, lines and panels are the same as in Fig¬ 
ure [2] We can see there is a general good agreement between the 
gravities and no significant offsets introduced by by the evolution¬ 
ary track method. 

nificant trends in the residuals were found as a function 
of reference mass, Teg, [Fe/H], V magnitude and tt (i?^ < 
0.12). Similar good agreements were also found for the 
radii and surface gravities. The analysis of the ages of 
binary systems provided an additional confirmation that 
the results obtained with the evolutionary tracks were 
consistent. 

The global results of our study suggest that the inter¬ 
polation of observed parameters in a grid of evolution¬ 
ary tracks, in particular the PARSEC set, is capable of 
providing accurate and relatively precise masses, radii 
and surface gravities for evolved stars with different ef¬ 
fective temperatures and metallicities. We acknowledge, 
however, that much work remains to be done both ob- 
servationally and theoretically in order to derive more 


Figure 6. Comparison between the ages of the members in binary 
systems. Symbols and lines are the same as in Figure[2] The upper 
panel shows a direct comparison of the ages. The lower panel shows 
the differences (in the sense A-B, relative to the mean age of the 
system and as a percentage) between the two sets of ages relative 
to the average age of the system. We can see there is a general 
good agreement between the ages. 

accurate input parameters, improve the models and bet¬ 
ter understand why some specific stars present a poor 
agreement in our comparisons. 

With this in mind, our team is currently working on a 
follow-up study that will repeat the analysis done here 
using different sets evolutionary tracks. This test could 
potentially determine the physics or free parameters that 
best reproduce the model-independent masses. We are 
simultaneously working on the detailed characterization 
of other potential benchmark stars using interferometry 
and asteros eismology, in a simila r fashion as was done for 
HD 185351 (j.Tohnson et al.l[2014f ). Finally, in the light of 
these new res ults, we plan to revisit the sample of retired 
A stars from [Johnson et ^ (|2010all . presenting new an¬ 
alyzes of their spectroscopic and kinematic properties. 

This research has made use of DEBCat and the SIM- 
BAD database, operated at CDS, Strasbourg, Erance. 
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Table 1 

Properties of the Benchmark Stars. 


Star 

V 

E{B - V) 

77 

Tea 

[Fe/H] 

t^max 

Au 




(mas) 

(K) 


(mHz) 

(mHz) 


Binary Stars 


AI Phe B 
CapellaA 
CapellaB 

9.542 ± 0.025 
0.892 ± 0.016 
0.763 ± 0.015 

0.020 ± 0.020 
0.000 
0.000 

5.78000 ± 0.37000 5010 ± 120 

75.99400 ± 0.08900 4970 ± 50 

75.99400 ± 0.08900 5730 ± 60 

-0.14 ± 0.10 
-0.04 ± 0.06 
-0.04 ± 0.06 



Stars with Asteroseismology 

11 Com 

4.740 ± 0.020 

0.016 ± 0.020 

11.25000 ± 0.22000 4841 ± 100 

-0.28 ± 0.10 

26.70 ± 1.34 

2.88 ± 0.06 

Arcturus 

-0.050 ± 0.010 

0.000 

88.83000 ± 0.54000 4286 ± 30 

-0.52 ± 0.04 

3.47 ± 0.17 

0.82 ± 0.02 

/3 Aql 

3.710 ± 0.009 

0.000 

73.00000 ± 0.20000 5030 ± 80 

-0.21 ± 0.07 

416.00 ± 20.80 

29.56 ± 0.10 

Note. — 

A portion of the 

table is shown 

here for guidance regarding its form 

and content. The complete table 

is available at 


https://drive.google.com/file/d/0B_C74xx43A0HTXBKMDhlYVI4RDA/view?usp= sharing 

AI P he B: Individual V magn itude was calcul ated from the value of in I Andersen et ay^^SSjj^consM ering t he reddening. E!(B — V) 

from fHrivnak & Milon^ | |1984n . Parallax from [Torres et al.| | |2010D . Tgff and [Fe/H] fro m[And ersen et al.| fTgSSl) . 

Capella A and B: Individual V magnitudes, parallaxes, TgffS and [Fe/H] from l^)rres et~^ ■I I I2015I1 . E(B — V) — 0.000 because d < 60 pc. 
11 Com: V m agnitude from SIMB AD. E{B — V) convert ed from Ay in |T^ed a et al.l 1120081) and a n arbitrary error of 0.02 was adopted. 
Parallax from f^n LeeuwenI <120071) . Teff and [Fe/H] from lTakeda et al.l l200^) ~ i^max and Ai/ from [Ando et all II2010D . Arbitrary errors 
of 5% and 2% in t^max and Ai.^, respectively. 

Arcturus: V magn itude from SIMBAD w ith an arbitrary error criV ) — 0. 01. E{B — V) — 0.0 00 because d < 60 p c. Parallax from 
Ivan LeeuwenI d2007l) . Teff and [Fe/H] from [Ramirez Sz Allende Prietd II2011D . i^max and Ai/ from fLagarde et al.l d2015l) . 

0 Aql: V magnitude from SIMBAD. E( B — V) — 0.00 0 because d < 60 pc. Parallax from fron LeeuwenI 120071) . Teff and [Fe/H] from 
IBruntt et al.l (120101) . i/'max and from fCorsaro et al.l (I2012D . Arbitrary error of 5% in i^max- 
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Table 2 

Evolutionary parameters for the benchmark stars. 


Reference Parameters PARAM Results 


Star 

M 

(Me) 

R 

(R@) 

logs 

M 

{M@) 

R 

(Re) 

logs 

Age 

(Gyr) 

References 


Binary Stars 

AI Phe B 
Capella A 
Capella B 

1.234 ± 0.004 
2.569 ± 0.007 
2.483 ± 0.007 

2.932 ± 0.048 
11.980 ± 0.570 
8.830 ± 0.330 

3.595 ± 0.014 
2.691 ± 0.041 
2.941 ± 0.032 

1.118 ± 0.104 
2.472 ± 0.104 
2.474 ± 0.036 

2.889 ± 0.246 
11.690 ± 0.285 
8.604 ± 0.212 

3.538 ± 0.089 
2.668 ± 0.027 
2.950 ± 0.021 

6.616 ± 2.122 
0.727 ± 0.117 
0.641 ± 0.017 

Torres et al. 

Torres et al. 

Torres et al. 

C2010~) 

(Ms) 

Stars with Asteroseismology 

11 Com 
Arcturus 
/3 Aql 

2.435 ± 0.423 
0.676 ± 0.121 
0.872 ± 0.133 

17.457 ± 1.141 
26.319 ± 1.829 
2.624 ± 0.134 

2.335 ± 0.022 
1.422 ± 0.022 
3.536 ± 0.022 

2.192 ± 0.330 
0.978 ± 0.057 
1.140 ± 0.105 

14.607 ± 1.212 
24.533 ± 0.663 
3.248 ± 0.125 

2.422 ± 0.128 
1.621 ± 0.040 
3,444 ± 0.066 

0.931 ± 0.381 
8.616 ± 1.653 
5.859 ± 1.868 

This work 
This work 
This work 



Note. — A portion of the table is shown here for guidance regarding its form and content. The complete table is available at 
https://drive.google.com/file/d/0B_C74xx43A0HTXBKNDhlYVI4RDA/view?usp=sharing 

^ The reference parameters are from the analysis of the RVs and light curves of the eclipsing binary system. The parameters derived from asteroseismology are: 
M — 1.717 zb 0.142 Mq , R — 11.480 zb 0.434 Rq and log g — 2.552 zb 0.006. The two sets of results agree within 1.6a. 
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Table 3 

Stars not Included in the Analysis. 


Star 


References 


Note 


ASAS 010538 B 
ASAS 182510 A,B 
ASAS 182525 A,B 
/3 Oph 
P UMi 
CF Tau A 
e Tau 
rj Her 
HD 50890 
HD 169751 
HD 170031 
HD 170053 
HD 175679 
L Dra 
Kepler-56 
Kepler-91 
KIC 3730953 
KIC 5737655 
KIC 6442183 
KIC 9705687 
KIC 11137075 
KIC 11674677 
M67 13 

OGLE SMC113.3 4007 A 
OGLE LMG CEP0227 B 
OGLE LMG-ECL-10567 A,B 


^^tajczal^^ah 

^atajcza^^^ah 

^atajcza^^^ah 

j4allinge^^^ 

Tarrant et al 


Ando et al. 

(20101 

Ando et al. 

(’ 20 T 0 I 





^racz^l^^aL 

Pileck^^al 
^ietrz^nsk^^ah 


[Fe/H] dependent on stellar evolution models 
[Fe/H] dependent on stellar evolution models 
[Fe/H] depends on stellar evolution models 
Error on reference mass >20% 

Error on reference mass >20% 

[Fe/H] dependent on stellar evolution models 
Suspicious (too high) reference mass 
Suspicious (too high) reference mass 
Error on reference mass >20% 

Error on reference mass >20% 

No independent parallax available 
Error on reference mass >20% 

Error on reference mass >20% 

No value for Ai/ 

No independent parallax available 
No independent parallax available 
Parallax error too large 
Error on reference mass >20% 

No independent parallax available 
Parallax error too large 
No independent parallax available 
Error on reference mass >20% 

Distance is not model-independent 
Assumed [Fe/H] 

[Fe/H] depends on pulsation models 
Remarks about the results for the system 


Note. — There are 11 K giants in IStello et al.l 1120081) . but the values of At' are not provided. 





































Masses of Retired A Stars 


15 


Table 4 

Statistics for the comparison between evolutionary track and reference parameters. 


Sample 

Stars 

(APar)^ 

{APar/Par He fT 


b'° 


(7 




(%) 






Par = M (M©) 


All Stars 

59 

-0.095 

± 

0.052 

-1.30 

± 

2.42 

0.60 

± 

1.65 

-3.76 

± 

4.79 

0.002 

1.32 

Binaries 

26 

0.041 

± 

0.065 

2.49 

± 

2.99 

0.57 

± 

1.93 

-2.24 

± 

5.76 

0.004 

1.34 

Asteroseismology 

33 

-0.202 

± 

0.073 

-4.29 

± 

3.58 

-11.53 

± 

2.93 

12.12 

± 

6.75 

0.333 

0.89 






Par = R 

{R&) 









All Stars 

59 

-1.002 

± 

0.172 

-4.81 

± 

1.32 

0.00 

± 

0.07 

-5.07 

± 

1.86 

0.000 

1.26 

Binaries 

26 

-1.033 

± 

0.214 

-2.57 

± 

0.73 

-0.11 

± 

0.04 

0.45 

± 

1.21 

0.268 

0.49 

Asteroseismology 

33 

-0.978 

± 

0.261 

-6.58 

± 

2.25 

-0.74 

± 

0.36 

-0.54 

± 

4.10 

0.122 

1.44 






Par = log 

g 









All Stars 

59 

0.010 

± 

0.011 

0.71 

± 

0.51 

-1.14 

± 

0.60 

3.04 

± 

1.77 

0.061 

1.22 

Binaries 

26 

0.003 

± 

0.015 

0.49 

± 

0.74 

-2.05 

± 

0.68 

4.74 

± 

1.70 

0.277 

0.75 

Asteroseismology 

33 

0.015 

± 

0.016 

0.89 

± 

0.70 

-1.28 

± 

1.14 

3.70 

± 

3.56 

0.040 

1.48 


^ APar — ParTrk. - , where Par can be mass, radius or surface gravity. The uncertainties are standard deviations 

of the means. 

^ Coefficients of the linear fit y = Ax + B. 












